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A system of two-dimensional weakly interacting bosons 
in a harmonic trap is considered. When the two-particle po- 
tential is a delta function Smith and Wilkin have analyti- 
cally proved that the elementary symmetric polynomials of 
particle coordinates measured from the center of mass are 
exact eigenstates. In this study, we point out that their 
proof works equally well for an arbitrary two-particle po- 
tential which possesses the translational and rotational sym- 
metries. We find that the interaction energy associated 
with the eigenstate with angular momentum L is equal to 
aN{N—l)/2+{b — a) NL/2, where a and b are the interaction 
energies of two bosons in the lowest-energy one-particle state 
with zero and one unit of angular momentum, respectively. 
Additionally, we study briefly the case of attractive quartic 
interactions. We prove rigorously that the lowest-energy state 
is the one in which all angular momentum is carried by the 
center of mass motion. 
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The study of Bose-Einstein condensation in trapped 
atomic gases has attracted a great deal of attention in the 
past few years . One of the central issues has been 
the possibility of creating quantized vortices in these di- 
lute atomic gases. Recently, vortex states in various sys- 
tems have been experimentally observed |||,^ . Some the- 
oretical investigations on quantized vortices and on ro- 
tating Bose condensates have also been carried out both 
in the Thomas-Fermi limit of strong interactions [7|-[lC|] 
and in the limit of weak interactions |11-18| between the 



atoms. For the weakly interacting bosons one is naturally 
led to study the model of N two-dimensional bosons in 
a harmonic trap with weak repulsive delta- function in- 
teractions Pt[| . One important theoretical problem here 
is to understand the properties of the system with a 
given total angular momentum L. Numerical studies 
have shown that when L > N the lowest-energy state 
is the one where an array of singly quantized vortices is 
formed |l5|,H. However, due to the complexity of these 
states, very few analytic results are known. In the range 
< L < N the structure of lowest-energy states turns 
out to be simpler. Numerical computations by Bertsch 
and Papenbrock showed that the interaction energy 
of the lowest-energy states has a very simple form and 
decreases linearly with L (the only exception being that 
for L = 1). They also noted that the wave functions 
for the lowest-energy states are simply the elementary 
symmetric polynomials of complex coordinates relative 
to the center of mass. Very recently, this remarkable 
formula for interaction ener gy was derived analytically 
by Jackson and Kavoulakis Moreover, Smith and 



Wilkin iQ proved that these symmetric polynomials are 
indeed exact eigenstates (see Ref. for an alternative 
proof). In this note, we point out that the analytic proof 
by Smith and Wilkin works not only for a delta-function 
potential but also for an arbitrary potential which pos- 
sesses translational and rotational symmetries. In this 
more general case, we find that the formula for the in- 
teraction energy still has the similar form of that for a 
delta-function potential. In addition, we briefly study 
the case of attractive quartic potential. We prove rig- 
orously that, as in the cases of attractive delta- function 
and harmonic potentials, the lowest-energy state is the 
one where all angular momentum is carried by the center 
of mass motion. 

Our model Hamiltonian is H = Hq + V , where 
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is the one-particle part and 
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is the two-particle interaction. Note that the two-particle 
potential has been taken to possess translational and ro- 
tational symmetries. We also assume that the interaction 
term is weak and does not make the system unstable. It is 
known that in the absence of interaction the one-particle 
spectrum is given by 



Enr.,1 = + 1^1 + 1, 



(3) 



where Ur is the radial quantum number and I is the an- 
gular momentum. For the system of N non-interacting 
bosons, the lowest-energy states with a given total angu- 
lar momentum L can be obtained by putting all bosons 
in the states with 71^ = 0, and with I of the same sign as 
L. Obviously, there is a huge degeneracy which equals 
to the number of ways to distribute L units of angular 
momentum among N bosons. 

When the two-particle interaction is added the degen- 
eracy will, in general, be lifted. If the interaction is 
weak enough, it suffices to analyze this problem by us- 
ing the first-order degenerate perturbation theory. We 
therefore need to diagonalize V in the restricted Hilbert 
space of lowest-energy states with angular momentum L. 
We shall do it in the second quantized form. For a pos- 
itive L (as we always assume), the relevant normalized 
one-particle wavefunctions for the states with = are 
given by 
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where / > and z = x + iy. is the complex coordinate. 
We denote by a; and aj the operators which annihilate 
and create one boson in the one-particle state Our 
Fock space is spanned by the basis states of the form 



(5) 



where |0) is the Fock vacuum and the occupation num- 
bers n'^s satisfy 



k k 
The second quantized form of V is 



aiak, 
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i,j,k,l 



where the matrix elements are given by 



Vijki — 



<fz<fz'r^ {Z) [Z') V{\Z - z'DV-fc (Z) {Z') . 



(8) 



Most of these matrix elements are actually vanishing. In- 
deed, the rotational symmetry of V implies that Vijki = 
unless i + i ~ k + I. This, of course, corresponds to the 
conservation of total angular momentum in two-particle 
collisions, an important feature of a rotationally symmet- 
ric potential. The translational symmetry of V also leads 
to a set of constraints on Vijki ■ Since the potential is real, 
this symmetry can be simply expressed by 



N 



d 



y—v 



= 0, 



(9) 



where = ^ — i-^^ ■ It is easy to check that the 
second quantized form of this condition reads 



(10) 



where = vAT+Tala/c+i. To work out the con- 



straints imposed by Eq.(|10D, we write 

V = Woojajaoflo + Via\a\aiai + vioajajaioo 

+Viia2aQaiai + vlia\a\a2ao + • • • • (11) 

Here we have intentionally displayed only five terms since 
the other terms turn out to be irrelevant for our analysis 
(see also Ref. @). Because V is hermitian, vo,vi and 
vio must be all real. For later uses, we also note that 



^'o — 2^0000; 



Vl = ^^1111- 



(12) 



Substituting Eq.(|rT|) into Eq.(|l^ yields 



0= L_,F 
= (uio - 2wo)aJaJaiao -I- {2vi 



Vio + V2vii)a\alaiai 
(13) 



Thus, the constraints, which are completely determined 
by the displayed terms in Eq.(pll), are 



(14) 



These two constraints are direct consequences of the 
translational symmetry and will play an crucial role in 
later analysis. 

To proceed further, it is necessary to introduce the 
elementary symmetric polynomials 



il<i2<---<iL 



(15) 



where L < N. The polynomial is related to the state 
\N~L,L) by 



{zi^Z2,- ■ ■ ,ZN \N - L,L) 



{N -L)\L\ 



1/2 
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Therefore, we shall use to represent the state 
\N — L, L). In fact, it is not hard to see that there 
is an one-to-one correspondence between the monomial 
symmetric polynomials of degree L and the basis states 
in our Fock space |]l8|. Consequently, each symmet- 
ric polynomial uniquely specifies a quantum state. We 
also like to remark that the coordinate representation 
of the operator = lI_ is simply NR, where R — 
'^j^Zi/N is the center of mass coordinate. Thus, the 
symmetric polynomial R^cl {k > 0) represents the state 
L'^\N — L, L). Furthermore, since commutes with 
V, we have vL'l \N - L, L) = L'^V \N - L,L). Equiva- 
lently, we may write it as 



VR''eL = R''VeL. 



(17) 



The elementary symmetric polynomials measured from 
the center of mass are defined by 



eM 



ll<l2<-<tL 



(18) 

Note that ei is the trivial zero function. It has been 
proven by Smith and Wilkin ||l^ that the state repre- 
sented by Eq.([l^) is an exact eigenstate when v{\r\) is a 
delta function. Here we follow the steps in their proof. 
First, we write 



M 



eM 



-EM) 



A/-L 



L=2 



{N - M)\{M - L)\ 



R^'-^CL 



^ ' {N-M)\M\ 



(19) 
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Next, by using Eq.(|ll|) we obtain 

V\N -L,L) 
= [vo {N- L){N- L-l)+ viL{L - 1) 
+vio{N - L)L]\N - L,L) 

+vii^L{L - 1){N - L + 1) |7V - L + 1, L - 2, 1) . (20) 

The monomial symmetric polynomial representing the 
state |iV - L + 1, L - 2, 1) is ig 



il<i2<---<iL-2 



NReL-i~LeL. (21) 



In terms of symmetric polynomials, Eq.(pO|) becomes 
VcL = [vo {N ~ L){N - L-l)+ viL{L - 1) 

+vw{N - L)L - -^viiL{N -L + l)]eL 
v2 

+ ^viiN{N-L + l)ReL-i. (22) 
V2 

It should be noted that even though Eq. ( pT| ) is valid only 
for L > 2 Eq.dl) holds in the whole range < L < N 
if the convention e_i = is taken. Up to this point, 
we have not made any use of the conditions given by 
Eqs.(|§. To uncover the role played by these conditions, 
we point out that Cm is a genuine eigenstate if Eq.(|2^) 
has the form: 

VeL = [fiN) + LgiN)]eL - [N - L + l)g{N)ReL-u 

(23) 

where / and g are two arbitrary functions of N . This 
statement can be easily proved by operating V on the 
right hand side of Eq.(p^) and by using Eqs.(p^) and 
( P3| ) . We also find that when Eq. (|2^) holds the eigenvalue 
associated with cm is given by 



SM^N^ f{N)+g{N)M. 



(24) 



It is remarkable that Eqs.(n4) are precisely the required 
conditions for making Eq.(|22|) of the form given by 
Eq.dl). Indeed, by substituting Eqs.(|l|) into Eq.(|2|) 
we get 



VeL = [voN{N - 1) + [vi ~ vo)NL\eL 
-{vi - vo)N{N -L + l)ReL-i. 



(25) 



As a result, bm is an exact eigenstate and the correspond- 
ing eigenvalue is 



eM,N = [voN{N - 1) + {vi - vo)NM] 

= \ [VomoN{N - 1) + (Fnii - Voom)NM] 



(26) 



Still we have a spectrum which varies linearly with the 
angular momentum M. It is worth noting that Vqooo 



(Villi) represents the interaction energy of two bosons 
in the lowest-energy state with I = Q (l = \). Eq.(|2^) 
represents the main result of the present work. As a ver- 
ification of this formula, we consider the delta-function 
potential 



V (r - 



27r?7(5 (r — r' 



(27) 



where 77 is a small dimensionless parameter. Simple cal- 
culations show 



1^0000 = m 



^1111 = ^'?- 



Consequently, Eq.(26) gives 



^delta 
-L,N 



N{N- 1) - -NL 



(28) 



(29) 



in agreement with the results of Refs. |13| , |l7| . 
Finally, we consider the case of quartic potential 

«(r-r') = ||r-r'|\ (30) 

When ?7 < the repulsive force between two particles 
which are moved away from the trap center in opposite 
directions would eventually be stronger than the trap- 
ping forces acting on them. We thus expect that in this 
case the harmonic trap is unable to stably confine the 
bosons in a finite region of space. To avoid instability 
we assume that the interaction is attractive;i.e., 77 > 0. 
The purpose of studying this case is two-fold. First, it 
serves as an additional example to which our results are 
applied. Secondly, it provides a new case, in addition 
to the cases of delta-function and harmonic potentials, 
where the lowest-energy states can be analytically deter- 
mined. After some straightforward algebras we get the 
following expression for the two-particle interaction: 

V = v[^N{N - 1) + ^-NL + h{L - 1) 

+ ^NJ-^{J+L^+L+J^) + h++L^^], (31) 

where N = J^k '^I'^fc; ^ — J2k ka\ak and we have defined 

k 

J+ = J1 = ^ (A: + 2) ^/k + lal^^ak] 

k 

L++ = lI_ = J2 v/(fc + 2)(fc + l)4+2afc- (32) 



In this case the interaction energy associated with cl is 
found to be 



quartic V 



N{N-1) + -NL 



(33) 
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The interaction energy increases linearly with L. Now cl 
is no longer the lowest-energy state with angular momen- 
tum L. We find that, as in the case of attractive delta- 
function potential , the lowest-energy state is the one 
where all angular momentum is absorbed by the center of 
mass motion. However, the proof used in Ref. |pl]] is not 
applicable here since the matrix elements Vijki are not all 
non-positive. We prove it by deriving a lower bound for 
the interaction energy. To this end, we first note that op- 
erating on an eigenstate amounts to adding an unit of 
angular momentum to the center of mass coordinate and 
the resulting state is still an eigenstate with the same 
eigenvalue. Therefore, only the states, called intrinsic 
states and denoted by |L, N)-^^, which have no center of 
mass excitation, need to be considered. Such states are 
characterized by 



0. 



(34) 



By using Eq.Q, {L, N\ L++L^^ \L, iV)„^, > 0, and 
„t {L, N\ J |L, 7V)^„, > L „u {L, N \L, iV),„, , (35) 
we obtain 



(L,iV|L,iV)„, 

1 



1 



^N{N -l) + NL+ -L{L - 1) 



(36) 



This inequality holds for all L > 0. Since |iV, 0) is the 
only eigenstate corresponding to eigenvalue r]N{N—l)/2, 
the lowest-energy state with a given angular momentum 
L must be 1-^,0), the state in which all angular mo- 
mentum is in the center of mass motion. It has been ar- 
gued ||ll| (see Ref. ||2^ , however) that this state is uncon- 
densed and is an example of the fragmented condensate of 
Nozieres and Saint James j2^. Moreover, the presence of 
such an uncondensed lowest-energy state is believed to be 
a general feature of the attractive interactions [0. Our 
study supplies one more example supporting this belief. 

In conclusion, we have extended the work by Smith and 
Wilkin to the models of N bosons in a 2D harmonic 
trap interacting via arbitrary rotationally and transla- 
tionally symmetric potentials. We have shown that bl 
remains an exact eigenstate and the associated interac- 
tion energy varies linearly with the angular momentum. 
Our analysis reveals the importance of rotational and 
translational symmetries in proving this result and makes 
clear the physical meaning of the coefficients appearing 
in the formula for the interaction energy. We also briefiy 
discuss the problem of attractive quartic interactions. A 
lower bound for the interaction energy of intrinsic states 
is derived. Based on this, we find that the lowest-energy 
state for a given angular momentum L is the one in which 
all angular momentum is carried by the center of mass 
motion. 

This work is supported by the National Science Coun- 
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Note added. A large class of interacting boson sys- 
tems is considered in a recent preprint by T. Papenbrock 
and G.F. Bertsch p^ . Our main results are in agreement 
with their findings. 
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